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ABSTRACT 

This  report  describes  the  improvements  made  to  the  ADINA  78  finite  element 
computer  code  in  order  to  permit  the  solution  of  structural  response  problems  involving 
very  rapid  heating.  A  material  model  was  added  to  the  two-  and  three-dimensional 
continuum  element  types  which  permits  the  modeling  of  temperature  dependent 
thermoelastic  anisotropic  material  properties.  A  different  material  model  was  also 
added  for  these  same  continuum  element  types  which  permits  the  modeling  of  strain  and 
temperature  dependent  thermoelastic  anisotropic  material  properties.  A  major  goal  of 
this  work  was  to  render  AD  IN  A  78  operational  on  the  ASC  computer  at  the  Naval 
Research  Laboratories  as  well  as  to  improve  the  computational  efficiency  of  the  code  by 
vectorizing  the  most  important  loops.  This  work  and  various  benchmark  test  cases  are 
also  described. 
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1.  INTRODUCTION 


For  the  past  several  years,  personnel  in  the  Materials  Science  and  Technology 
Division  of  the  Naval  Research  Laboratories  (NKI.)  in  Washington,  U.  C.  have  been 
actively  engaged  in  research  directed  at  understanding  and  predicting  the  effects  of 
rapid  heating  of  naval  structural  components.  These  efforts  quickly  focused  on  a  small 
number  of  general  purpose  finite  element  computer  programs  that  were  capable  of 
solving  these  problems.  The  ADINA^1'^  code  was  selected  for  these  studies  because  it 
was  the  only  general  purpose  code  for  which  source  coding  was  readily  available  at 
reasonable  cost,  and  because  the  source  coding  was  well  written,  modular,  well 
documented  and  easy  to  modify  for  special  applications.  This  report  describes  the 
improvements  made  to  ADINA  78^1  in  order  to  make  it  more  suitable  for  the  solution  of 
structural  response  problems  involving  rapid  heating. 

The  1978  version  of  ADINA  contained  a  material  model  (Model  3)  for  isotropic 
thermoelastic  temperature  dependent  properties  and  a  material  model  (Model  2)  for 
anisotropic  elastic  properties  for  the  two-  and  three-dimensional  (2D  &  3D)  solid 
elements.  However,  there  was  no  anisotropic  thermoelastic  temperature  dependent 
model.  This  capability  was  added  to  ADINA  78  using  Model  12  for  the  2D  &  3D  solid 
element  types.  The  implementation  of  this  model,  the  necessary  data  cards  together 
with  sample  problems  are  described  in  Section  2  of  this  report. 

The  main  thermal  problems  of  interest  to  NR  L  involve  temperature  excursions  of 
hundreds  of  degrees  (°F),  consequently,  the  material  response  becomes  nonlinear.  h 
order  to  approximately  model  these  thermal  nonlinearities  tor  anisotropic  materials. 
Model  6  was  also  added  to  ADINA  78  for  the  2  D  &  3D  solid  types.  The  implementation  of 
this  material  model,  the  necessary  data  cards  together  with  sample  problems  are 
described  in  Section  3  of  this  report. 

Recently,  the  Naval  Research  Laboratory  acquired  a  new  and  unique  fourth 
generation,  high  speed  digital  computer,  a  Texas  hstruments  Advanced  Scientific 
Computer  (  ASC).  While  the  ASC  is  no  longer  in  production  and  although  there  are  only  a 
few  left  in  operation,  it  is  a  very  powerful  computer.  The  rated  speed  of  the  ASC  for 
scalar  operations  is  approximately  equal  to  a  CDC  6800  computer;  additionally,  the  ASC 
has  special  hardware  and  software  capabilities  to  perform  ‘vector’  floating  point 
operations  substantially  faster  than  scalar  operations.  Prior  to  the  current  work,  ADINA 


» 
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had  been  run  only  in  the  scalar  inode.  Thus,  a  further  facet  of  the  current  work  was  to 
"vectorize*  ADINA  to  improve  its  computational  efficiency  on  the  ASC.  This  work  and 
various  benchmark  test  cases  are  described  in  Section  4. 


Section  5  gives  brief  conclusions  regarding  the  current  effort  and  Section  6 
contains  recommendations  for  further  improvements  and  modifications  of  ADINA. 

Three  appendices  have  been  included  to  complete  the  report.  Appendix  A 
describes  the  ADINA  files  on  the  ASC,  Appendix  B  briefly  discusses  various  execution 
problems  with  ADINA,  and  Appendix  C  describes  an  /Anisotropic  Deformation  Theory 
Plasticity  that  was  considered  for  ADINA. 


2.  TEMPER  ATURE  DEPENDENT  AN  ISO  1  R  OP  1C  I  HE  RMOE  LAST  1C  MODE  L 


2.1  General 


I  he  70  version  of  ADINA  lias  available  a  large  number  of  material  models  for  2  D 
and  313  solid  element  types  as  indicated  in  Table  2.1. 

Table  2.1  ADINA  78  Material  Models  for  2D  &  3D  Solid  Element  Types 


NP  AR(1 5) 
Model  No. 
1 
2 

3 

4 

5 

7 

8 

9 

10 
11 
13 


Model  Name 

Isotropic  Linear  Elastic 

Orthotropic  (  Anisotropic)  Linear  Elastic 

Isotropic  Thermoelastic 

Curve  Description 

Concrete 

Elastic-Plastic  (  Drucker-P  rage r-  Cap) 

Elastic-Plastic  (von  Mises,  Isotropic  Hardening) 

Elastic-Plastic  (Kinematic  Hardening) 

Thermoelastic-Plastic-Creep  (von  Mises,  botropic  Hardening) 
Thermoelastic-Plastic-Creep  (von  Mises,  Kinematic  Hardening) 
hcompressible  Nonlinear  Elastic  (Mooney-Rivlin),  2  D  elements  only 


Unfortunately,  ever  since  the  original  version  of  ADINA,^'^  bathe  has  retained  this 
cumbersome  identification  of  the  material  models.  For  example,  while  ADINA  can 
model  both  orthotropic  and  thermoelastic  problems  via  Model  2  and  3,  respectively,  it 
does  not  permit  the  modeling  of  orthotropic  thermoelastic  problems  although  the 
necessary  coding  for  the  calculation  of  the  thermal  force  vector  and  local  axes  of 
material  orthotrophy  exist. 

In  order  to  properly  model  typical  naval  structures  subjected  to  rapid,  intense 
heating,  it  is  necessary  to  have  an  anisotropic  thermoelastic  temperature  dependent 
material  model.  Under  the  current  effort  this  capability  was  added  to  ADINA  78  by 
implementing  a  new  Model  12  as  described  in  the  remainder  of  this  section. 
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In  terms  of  compatibility  with  other  ADINA  options.  Model  12  is  identical  to 
Model  3.  It  will  work  for  dynamic  or  static,  geometrically  linear  or  nonlinear  analyses, 
etc.  However,  Model  12  is  intended  for  static,  materially  nonlinear  analyses  only.  Its 
use  with  other  options,  while  allowable  as  input  to  ADINA,  may  produce  erroneous 
results.  The  user  is  encouraged  to  carefully  verify  use  of  Model  12  with  these  options. 
This  was  not  done  under  the  current  effort. 

Model  12  is  very  similar  to  Model  3  in  the  method  used  to  store  and  retrieve 
material  properties  as  a  function  of  temperature  from  the  PROP  array  (see  Section 
2.3).  The  strain  and  temperature  dependent  Model  6  described  in  Section  3  is  likewise 
similar  to  Models  3  and  12.  While  Model  b  can  be  used  without  strain  dependency  to 
model  temperature  dependent  anisotropic  properties,  this  usage  is  not  recommended; 
Model  12  is  intended  for  these  problems  and  is  easier  to  use  and  more  efficient. 

2.2  Implementation  of  MO DE  L  12 

Because  ADINA  was  capable  of  treating  isotropic  thermoelastic  materials  and 
anisotropic  materials,  the  implementation  of  Model  12  was  straightforward,  albeit 
nontrival.  This  section  briefly  describes  the  modifications  and  additions  that  were  made 
to  ADINA. 

Logic  was  added  in  the  element  calling  routine  for  the  2D  (TODMFE)  and  3D 
(THRtDM)  solid  element  types  to  establish  storage  boundaries  for  the  new  material 
property  information  that  must  be  stored  for  Model  12.  Logic  was  modified  in  the 
M  AT  KT2  (2D)  and  M  ATW  RT  (3D)  subroutines  to  read  and  print  the  new  material  property 
inf  nation.  Minor  changes  were  made  to  the  2D  routines  TDFE  and  QUADS  in  order  to 
recalculate  the  orientation  of  the  local  material  axes;  no  similar  changes  were  needed 
for  the  3D  elements. 

h  the  2D  and  3D  solid  element  overlays  there  are  subroutines  present  for  each 
material  model  (1-13).  Dummy  subroutines  EL2D12and  EL3D12for  Modell  2  existed  and 
these  were  replaced  with  functional  logic  that  called  the  following  suite  of  new 
subroutines  that  were  added: 


2D 


3D 


ITH212 
I  11212 

M KO  I  2D 


IT  1131 2 
T 113 12 
M  KOT3D 


While  the  corrections  were  straightforward,  they  were  nontrivial  and  fairly  extensive  in 
that  approximately  1000  lines  of  code  were  added  or  modified. 

2.3  Data  Cards  for  MODE  L  12 

A  small  number  of  modifications  to  the  ADINA  7ii  User's  Manual^  are  necessary 
in  order  to  describe  how  to  use  Model  12.  However,  there  are  no  changes  for  any  other 
material  models.  Rather  than  reproducing  extensive  portions  of  the  user’s  manual,  in  the 
sequel  only  the  modifications  and  additions  to  the  manual  are  described. 

The  input  requirements  for  Model  12  are  essentially  the  same  as  for  the 
hypothetical  union  of  Models  2  and  3  and  they  follow  the  same  format;  they  could  be 
deduced  from  the  existing  manual  by  a  clever  user.  The  one  difference  is  actually  an 
improvement  designed  to  eliminate  needless  input;  M=NPAK(14)  on  the  Element  Croup 
Control  card  has  been  used  to  define  the  number  of  temperatures  at  which  the  user  wants 
to  input  different  properties.  At  least  2  discrete  temperatures  should  always  be  given 
(M  >  2).  There  is  no  practical  upper  I imit  on  M. 

The  material  properties  are  stored  in  a  table  in  a  manner  similar  to  the  three- 
dimensional  array  PROP  described  on  the  succeeding  pages.  During  the  solution  phase, 
the  properties  are  linearly  interpolated  from  the  table  as  needed  as  illustrated  in  Figure 
2-1.  The  temperatures  must  be  given  in  strict,  discrete,  ascending  order  (T-j  <  T2  <  ...  < 
T  );  no  two  can  be  the  same.  Note  that  the  properties  are  taken  as  undefined  outside 
the  given  temperature  range  and  execution  is  terminated  if  T-j  >  T  >  T ni»  The  strain- 
stress  law  and  its  inverse  the  stress-strain  compliance  law  are  shown  in  Figures  2-2  and 
2-3,  respectively.  For  Model  12  the  interpolation  is  performed  on  the  stress-strain 
compliances  given  in  Figure  2-3  and  not  on  the  moduli,  Poisson's  ratio  and  coefficients  of 
thermal  expansion. 

In  order  to  avoid  needless  duplication  within  this  report,  the  following  description 
combines  the  data  cards  for  the  2D  (Section  X  1 1  of  the  AD  IN  A  7b  Use  r's  Manual^  pages 
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Temperature 

Model  12  (and  6)  Temperature  Interpolation  Method. 
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where  the  Gj  and  Oj  are  tile  normal  strains  and  stresses  in  t he  Xj  directions,  and  y j j  and  Ojj 
are  tlie  corresponding  engineering  shear  strains  and  shear  stresses, respectively,  and  A  T  is 
the  temperature  above  the  stress  free  refrence  temperature,  A  1  =  T-T  R  E  F. 


Note  v.  /E  -  v  ,/E,,  etc. 
b  a  a  a  b  b 


Figure  2-2  Strain-Stress  Law  for  MODELS  3  and  12. 
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Figure  2-3  Stress-Strain  Compliances  for  MODE  LS  3  and  12. 
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4  -  17)  ami  3D  (Section  X  lllof  the  manual,  pac.es  4  -  17)  solid  element  types. 
Modifications  for  the  2D  and  3D  Solid  Element  Input  Data  (MODE  L  12) 


Columns 


Variable 


Element  Croup  Control  Card  (2014) 
Entry 


53-56 

57-60 


M=NP  AK(  14)  Number  of  temperature  points  (Models  6&12  only) 

NP  AR(15)=12  Material  model  number;  12  is  the  Orthotropic 

I  hermoelastic  Model  number 


Set  NPAR(15)  =  12  to  activate  material  Model  12  input  data.  M  =  NP  AR(1 4)> 2 
must  also  be  set  for  material  Model  12.  Otherwise,  the  Element  Group  Control  Card 
data  is  unchanged. 

Material  MODEL  12  Property  Data  Cards 

Each  of  the  following  card  sets  must  have  at  least  1  card  with  2  temperature 
points  (M  >  2).  If  M  >  8,  use  as  many  cards  as  needed  in  format  8 FI 0.0.  The  2D  and  3D 
solid  material  model  12  requires,  respectively,  10  and  12  card  sets  corresponding  to  the 
following  properties:  Ea;  Efa;  Ec;  vab;  vac;  vbc;  Cab;  Cac*;  Cbc*;  «a;  «b;  and  v  The 
subscripts  a,b,c  refer  to  the  local  orthogonal  principal  material  axes  as  defined  in  Section 
X  II  and  XIII,  respectively,  for  2D  and  3D  solids.  For  2  D  solids,  the  a-b  axes  lie  in  the  y-z 
plane  of  modeling.  Note  that  while  the  stress-strain  law  is  orthotropic  in  the  local  (a,b,c) 
material  coordinate  system  (Figure  2-3),  it  is  “anisotropic’  in  the  global  coordinate 


*  No  t«  that  G  and  Gbc  arc  not  defined  for  2D  solids;  thus,  card 
sets  8  and  9  are  om itted  for  2D  solid  el  erne  n  t  s . 
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system.  A  local  anisotropic  model  could  easily  he  added  at  the  expense  of  4  additional 
data  sets  for  2D  and  15  additional  sets  for  3D  solids.  The  following  describes  the 
material  property  data  cards  read  for  Model  12. 


Columns  Variable 


Card  Set  1 

I- 10  PROP(1,N,1) 

II- 20  PROP(2,N,1) 


•  • 

PK0P(M,  N,1) 

Card  Set  2 

I- 10  PK0P(1,N,2) 

II- 20  P  K  OP(2,  N,2) 


•  • 

PR0P(M,N,2) 

Card  Set  3 

I- 10  PKOP(1,N,3) 

II- 20  PROP(2,N,3) 


•  » 

PKOP(M,N,3) 

Card  Set  4 

I- 10  PROP(1,N,4) 

II- 20  PROP(2,N,4) 


P  R  0P(M ,  N,4) 

Card  Set  S 

1-10  PK0P(1,N,5) 

11-20  P  K  Ol’(2,  N,5) 


PROP(M,N,5) 


Lnt  ry 

1C,  temperature  at  point  1 
T 2  temperature  at  point  2 

• 

1  m  temperature  at  point  m 

ta1  a'direction  modulus  at  point  1 
^a2  a_direction  modulus  at  point  2 

• 

^  am  a“dircction  modulus  at  point  rn 

Eb1  b"c,i  rection  modulus  at  point  1 
Cb2  ^-direction  modulus  at  point  2 

• 

Ebm  l>_di rection  modulus  at  point  m 

Ec1  direction  modulus  at  point  1 
Ec2  c-c*'  rection  modulus  at  point  2 

• 

Ecm  c"direction  modulus  at  point  m 

vab1  ai>'(>°'sson's  ratio  at  point  1 
Vab2  ab"*,o<sson's  ratio  at  point  2 

• 

vabm  al,“Po<sson's  ratio  at  point  m 
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Card  Set  0 

1-10 

11-20 

• 

PKOP(1,N,6) 

PKOP(2,N,6) 

uac1  ac-1>°'sst,n's  ratio  at  point  1 
vac2  ac“l>u'sson's  ratio  at  point  2 

• 

« 

I’K  OP(M,N,6) 

• 

v  ac-Poisson's  ratio  at  point  in 

Card  Set  7 

1  -  10 
11-20 

• 

PROP(1,N,7) 

1'  l\  OP(2,  N,7) 

• 

vbr1  Dc-Poisson's  ratio  at  point  1 
vbc2  De-Poisson's  ratio  at  point  2 

« 

• 

• 

• 

• 

PROP(M,N,7) 

• 

vbcm  De-Poisson's  ratio  at  point  rn 

Card  Set  0* 

1  -  10 
11-20 

• 

PROP(1,N,«) 

PROP(2,N,«) 

• 

Gabi  ab-shear  modulus  at  point  1 

Cab2  ab-shear  modulus  at  point  2 

• 

• 

« 

• 

• 

PROP(M,N,8) 

• 

• 

Ci  k  ab-shear  modulus  at  point  rn 

Card  Set  9* 

1-10 

11-20 

• 

P  ROP(1,N,9) 

P  K  OP(2,  N,9) 

• 

G  .  ac-shear  modulus  at  point  1 
t»aC2  ac-shear  modulus  at  point  2 

• 

• 

• 

• 

% 

PROP(M,N,9) 

• 

• 

G  m  ac-shear  modulus  at  point  in 

Card  Set  10 

1  -  10 
11-20 

• 

PP.  OP(1,N,10) 

P  i<  OP(2,  19,1 0) 

• 

^bcl  De-shear  modulus  at  point  1 

C^c2  De-shear  modulus  at  point  2 

• 

• 

• 

• 

• 

PROP(M, 19,10) 

• 

• 

^bem  De-shea r  modulus  at  point  m 

Card  Set  11 

1  -  10 
11-20 

• 

P  ROP(1,N,11) 

P  K  OP(2,N,1 1) 

• 

a-direction  expansion  coefficient  at  point  1 
aa2  a"^irection  expansion  coefficient  at  point  2 

9 

• 

• 

• 

• 

PROP(M,N,11) 

• 

• 

aama-direction  expansion  coefficient  at  point  m 

•Card  Sets  0  ft  9  must  be  omitted  for  two-dimensional  problems. 


»  a 


ID 


PK0P(1,N,12) 
P  R  OP(2,  N,1 2) 


ct^i  i)-direction  expansion  coefficient  at  point  1 
ftb2  biii rection  expansion  coefficient  at  point  2 


Card  Set  12 

1  -  10 


11-20 


PKOP(M,N,12)  a^b-di rection  expansion  coefficient  at  point  m 

Card  Set  13 

1-10  PUOP(1,N,13)  «ci  c-direction  expansion  coefficient  at  point  1 

11-20  PKOP(2,N,13)  £tC2  c-direction  expansion  coefficient  at  point  2 


Card  Set  14 

1  -  1u 


PR0P(M,N,13)  acmc-direction  expansion  coefficient  at  point  m 
PKOP(1,N,14)  TREF  reference  stress  free  temperature 


».  v 


2.4  Sample  Problems  for  MODEL  12 

A  two-dimensional  and  similar  three-dimensional  sample  problem  was  run  to 
validate  the  new  Model  12  logic  added  to  ADIiMA  7b.  Figures  2-4  and  2-5  show  the  2D 
atid  41)  problem  geometries,  respectively.  The  20  problem  is  a  plane  strain  90°  segment 
of  a  thin  axisymmetric  ring.  The  mean  radius  of  the  ring  is  20'  and  the  thickness  is 

0.5*. 

Symmetry  boundary  conditions  are  applied  on  the  y-(uz=0)  and  z-(u^=0)  axes,  so 
that  the  solution  is  symmetric  about  the  x-axis.  A  total  of  nine  8-node,  bi-quadratic 
quadrilateral  elements  spaced  at  10°  and  4G  nodal  points  were  used  as  shown  in  Figure 
2-4.  The  loading  consisted  of  a  333  psi  external  pressure  and  the  thermal  loading 
described  below. 

Temperature  dependent  orthotropic  moduli  and  orthotropic  thermal  expansion 
coefficients  were  used  as  shown  in  Figures  2-6a  &  b.  The  stress  free  reference 
temperature  was  TREF  =0  and  the  temperature  at  which  the  calculation  was  done  was 
T  =  10.  Thus,  from  Figures  2-6a  &  b 


E  =  lx  10°  psi 

d 

Eb  =  1.5  x  106  psi 
Ec  =  2  x  106  psi 


Gab=  -1 

x  10 

psi 

«a  =  .01 

t 

aaAT  = 

%  =  '°2 

t 

%AT  = 

ac  =  *03 

t 

«CAT  = 

.1 

.2 

.3 


Poisson's  ratios  were  taken  to  be  temperature  independent,  however,  they  were 
orthotropic 

vab  =  -1 
vac  =  *2 
vbc  =  '3 

The  ADINA  78  data  cards  are  given  in  Tables  2.2  and  2.3  for  the  2D  and  3D 
problems,  respectively.  The  2D  Model  12  problem  corresponds  to  the  file  R2D3  in 
Appendix  A  and  the  3D  problem  to  file  K3D2.  The  solutions  to  the  2D  and  3D  problems 
arc  identical  (in  closed  form,  not  numerically)  and  are  as  follows 
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333  psi  External  Pressure 

/ 

\/ 


2  2 
1  2 


0. 

0. 

10. 

.1  0  . 

2  2 

0. 

10. 

10. 

10. 

X 

1  1 

0 

1 

1 

1 

1 

X 

4  1 

0 

0 

1 

1 

1 

X 

44  1 

0 

0 

1 

1 

1 

X 

6  .1 

0 

0 

1 

1 

1 

X 

46  1 

1 

0 

1 

1 

1 

X 

2  1 

0 

1 

1 

1 

1 

X 

7  1 

0 

0 

1 

1 

1 

X 

47  1 

1 

0 

1 

1 

1 

X 

3  1 

0 

1 

1 

1 

1 

X 

5  1 

0 

0 

1 

1 

1 

X 

45  1 

0 

0 

1 

1 

1 

X 

0  1 

0 

0 

1 

1 

1 

X 

40  1 

1 

0 

1 

1 

1 

1 

1 

• 

1 

10. 

1 

40 

10. 

2  9 

1 

1 

1 

2. 

20. 

2. 0E5 

2 . 

026 

3. 0E5 

3. 

BEG 

4. 0E5 

4  . 

BEG 

.  1 

.  1 

.  2 

.  2 

.  3 

.  3 

2. 0E4 

2. 

0E5 

.  002 

.02 

.  004 

.  04 

.  006 

.  0S 

0. 

1  B 

1 

1 

5 

, 

007266463 

3  1 

6 

8 

2 

4 

7  5 

9  B 

1 

1 

5 

087266463 

43  41 

46 

43 

42 

44 

47  45 

2  1 

1 

6 

4 

32 

3.333333 

333.33 

333 

1 . 0 

2  1 

41 

4G 

44 

*:» 

3.  333333 

333.33333 

1  . 

1  2 

1  . 

1 

40  2 

1  . 

£>  TOP 


20 . 25 

20 . 25 

5. 

20.25 

05. 

20.25 

10. 

20.25 

90. 

20.00 

20.00 

IB. 

20.00 

98. 

19.75 

19.75 

5. 

19.75 

B5 . 

19.75 

10. 

19.75 

90. 

2 

12 

i  .  0 
1 . 0 


Table  2.2.  2D  Solid  Element  Model  12  Sample  Problem  -  File  R2D3 
Appendix  A. 


97100m 

2 


-  3D  PLANE  STRAIN  MODEL  12 

111  1  . 


10.  ?5 

19.75 

5 . 

10 

19.  75 

B5 . 

19.75 

10. 

10 

19. 75 

90. 

20 . 00 

20.00 

10. 

10 

20.00 

90. 

20.25 

20.25 

5 . 

10 

20.25 

05. 

20.25 

10. 

10 

20 . 25 

90. 

-1 . 0 

19.75 

“1 . 0 

19.75 

5. 

10 

-1.0 

19.75 

05. 

-1 . 0 

19.  75 

10. 

10 

-1.0 

19.75 

90. 

-1.0 

20.00 

-1.0 

20.00 

10. 

10 

-1 . 0 

20. 00 

90. 

-1 . 0 

20.25 

-1.0 

20.25 

5. 

10 

-1.0 

20.25 

85. 

-1.0 

20.25 

10. 

10 

-1.0 

20.25 

90. 

2  12 


Table  2.3.  3D  Solid  Element  Model  12  Sample 
Problem  -  File  R3D2  Appendix  A. 
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1 

97 

7 

95 

1 

97 

17 

95 

1 

97 

27 

95 

1 

97 

37 

95 

1 

97 

47 

95 

i 

97 

57 

95 

1 

97 

67 

95 

l 

97 

7  7 

95 

1 

97 

87 

95 

1 

16 

16 

1 

1 

1 

10 

1  o 

15 

5 

6 

12 

1 J 

l 

2 

8 

9 

0 

10 

0 

7 

0 

r  i 
(.• 

14 

13 

4 

0 

1. 

16 

16 

1 

1 

2 

10 

vi  r. 

25 

15 

16 

22 

21 

1 1 

12 

n 

19 

0 

20 

G 

17 

0 

J  0 

24 

23 

13 

14 

0 

3 

16 

16 

1 

1 

3 

10 

36 

35 

25 

26 

32 

31 

21 

22 

0 

39 

0 

30 

0 

27 

0 

28 

34 

33 

23 

24 

O 

4 

16 

16 

1 

1 

4 

10 

46 

45 

35 

36 

42 

41 

31 

32 

0 

39 

0 

40 

0 

37 

0 

38 

44 

43 

33 

34 

G 

5 

16 

16 

1 

1 

5 

10 

56 

55 

4  5 

46 

52 

51 

41 

42 

0 

A  49 

0 

50 

0 

47 

0 

40 

54 

53 

43 

4  4 

0 

6 

16 

*  16 

1 

1 

6 

10 

66 

65 

55 

56 

62 

61 

51 

52 

0 

59 

0 

60 

0 

57 

0 

58 

64 

63 

53 

54 

0 

7 

16 

16 

1 

1 

7 

10 

76 

75 

65 

66 

72 

71 

61 

62 

0 

69 

0 

70 

0 

67 

0 

68 

74 

73 

63 

64 

B 

0 

16 

16 

1 

1 

8 

10 

86 

85 

75 

76 

82 

81 

71 

72 

0 

79 

0 

80 

0 

77 

0 

78 

04 

83 

73 

74 

0 

9 

16 

16 

1 

1 

9 

10 

96 

95 

85 

86 

92 

91 

81 

82 

0 

89 

0 

90 

0 

87 

0 

80 

94 

93 

83 

84 

D 

1 

16 

15 

5 

6 

0 

9 

0 

10 

333 .  7:333 

333.3333 

333.3 

333  333.33 

33 

10 

1 

96 

95 

85 

86 

0 

89 

0 

90 

333.3333 

333.333 

3  : 

333.3333  333.33 

33 

1 

2 

1 . 

1 

97 

2 

1 . 

STOP 

Table 

2.3. 

(cont'd) 

3D 

Sol  id 

Element 

Model 

12  Sample  Problem 

-  File 

R3D2  Appendix 
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where  r,  the  radial  direction,  corresponds  to  the  AUINA  a  direction,  x  corresponds  to  the 
c  direction  and  b  to  the  U  or  hoop  direction.  Note  the  2  D  properties  have  been  rotated  5° 
(.07.1  radians)  so  that  a  is  in  the  radial  direction  at  the  center  of  the  element.  (The  I1L1 
material  angle  listed  on  page  XI  1.43  of  the  ADINA  manual  must  be  given  in  radians,  not 
degrees  as  the  manual  specifies.)  Likewise,  the  3D  material  axes  orientation  cards 
rotate  the  properties  so  that  the  a  direction  is  in  the  radial  direction  at  the  center  of  the 
element.  Thus,  the  results  are  numerically  symmetric  about  the  x-axis.  The  strain- 
stress  law  becomes 


because  the  ring  is  thin  and  free  to  expand,  the  hoop  stress  is 

a  -  -  =  -13333  psi 

0  t 

Solving  for  the  3  remaining  unknowns  gives 

o  =  -606000  psi 

x 

and 

e  r  =  .1615 
eG  =  .2820 

Thus  the  radial  displacement  is 

u  =  It  c,  =  5.640 

I  J 

The  numerical  solutions  given  by  ADINA  agree  closely  with  these  values  as  shown 
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in  Tables  2.4  and  2.5 


T 


1  able  2.4 

Computed  Nodal 

Point  Displacements 

2D 

3D 

Kadius 

Node 

Displacement 

Node 

Displacement 

19.75 

3 

5.597 

1 

5.597 

20 

2 

5.639 

3 

5.639 

20.25 

1 

5.680 

5 

5.680 

Table  2.5  Computed  Element  Stresses 


2D 

3D 

Output  Pt. 

2  D/3  D 

°x 

°x 

°0 

1/2 

-604504 

-14581 

-604504 

-14835 

3/3 

-603653 

-11917 

-603653 

-12128 
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3.  TEMPER  AIUKE  DEPENDENT  AN  ISO  X  KOP  1C  TH  E  R  M  OE  L  AST  1C  MODE  L 


3.1  General 

Over  a  temperature  range  of  several  hundreds  of  decrees  (°  I  ),  the  materials  used 
in  important  naval  structures  experience  a  significant  level  of  strain.  The  mechanical 
properties  of  these  same  materials  are  altered  .as  a  (unction  of  the  level  of  strain  as  well 
as  a  function  of  the  temperature.  The  variation  in  mechanical  properties  as  a  function  of 
temperature  can  be  modeled  in  ADINA  7iS  as  described  in  Section  2.  However,  there  is 
still  a  need  to  model  the  variation  in  properties  as  a  function  of  strain.  This  section 
describes  the  addition  to  ADINA  7U  of  material  Model  6  that  is  capable  of  modeling  the 
strain  dependency  and  the  combined  strain  and  temperature  dependency  of  the  material 
properties. 

Physically  Nonlinear  Elastic  Materials 

The  dependency  of  the  elastic  constants  on  strain  is  a  class  of  problems  in 
mechanics  known  as  physically  nonlinear  or  nonlinear  elasticity  problems.  This  class 
distinction  is  independent  of  kinematic  nonlinearities  (large  deformations  or  finite 
elasticity)  and  temperature  dependency.  It  is  an  extremely  c  inplex  class  of  problems 
that  is  not  frequently  found  in  the  literature^'5^  and  often  is  improperly  treated 
(especially  in  soil  mechanics).  Tire  proper  invariant  treatment  of  materials  of  this  class 
requires  the  development  of  a  constitutive  equation  of  the  form 


where  U  =  U  ( I-,,  ^  *3.  A^)  is  the  strain  energy  density  expressed  as  a  function  of  the  3 
invariants  of  the  strain  tensor 


cii 


Eiieji 


ijcjkf  ki 
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ami  the  3  vectors  (  Aj,)  describin'*  the  local  principal  material  directions  for  anisotropic 
materials.  Note  that  even  if  the  strain  energy  density  does  not  depend  on  A^,  the 
material  is  isotropic  only  about  the  undeformed  state  (  L  =  I  =  I  ==0)  and,  in  general,  is 
highly  anisotropic  otherwise. 

I  he  strain  energy  density  for  isotropic,  linear  elastic  materials  is  expressed 


where  p  and  A  are  the  Lame  constants.  Thus 


u  .  .  =  2 1 1  c  .  . 

•I  i  I 


A  e 


kk 


6  . 

i  J 


For  a  typical  nonlinear  elastic  material,  tests  must  be  performed  to  measure  the 
unknown  "elastic”  constants  in  a  polynomina!  expansion  of  the  strain  energy  density 


U  - 


A  I  \  * 
12  1 
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*  2  + 


1  3 


A2  3 
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*3  + 


A 


14  1 


2  4  1 


I,  I. 


A 


3  4  1 


3  S 
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for  even  the  simplicst  material  this  task  becomes  hopeless.  There  are  so  many  constants 
that  an  almost  infinite  number  of  tests  must  be  performed  to  determine  their  values. 
This  type  of  approach  has  rarely  been  successful.  (However,  one  exception  is  the 
Mooncy-Rivlin  rubber  like  material  model  available  in  ADINA.) 

A  more  tractable  approach  is  the  deformation  theory  plasticity  model  described  in 
Appendix  C.^l  Deformation  theory  plasticity  seeks  to  approximate  nonlinear,  path 
independent,  stress-strain  behavior  for  loading  paths  that  experience  no  unloading.  These 
theories  have  been  successfully  used  to  treat  several  anisotropic  materials  commonly 
used  in  naval  structures. Deformation  theory  requires  the  solution  of  only  a  path 
independent  nonlinear  elasticity  problem,  not  a  plasticity  problem.  ADINA  has  a 
nonlinear  static  solution  procedure  that  uses  a  secant  stiffness  matrix  that  is  very 
compatible  with  the  deformation  theory  approach. 

Implementation  of  the  deformation  theory  plasticity  model  described  in  Appendix 
C  in  the  ADINA  code  would  have  required  a  major  effort  well  beyond  the  scope  of  the 
present  work.  Such  an  effort  would  be  justified  only  if  the  theory  were  essential  to 


2? 


predict  iiij;  the  response  of  structural  materials  to  rapid,  intense  treating.  It  was  felt  that 
this  level  of  effort  v\.is  not  justified,  and  instead  it  was  decided  to  implement  asimplier 
model  that  would  permit  the  user  to  describe  the  usual  anisotropic  elastic  constants  as  a 
function  of  strain.  While  this  model  is  not  invariant,  it  can  lx;  used  to  approximately 
model  mild  nonlinearities  in  materials  that  undergo  path  independent  loading  only.  I  he 
remainder  of  this  section  describes  the  material  Model  0  that  was  implemented  in  ADINA 
to  treat  strain  and  temperature  dependent  properties. 

3.2  Implementation  of  M ODt  L  fa 

The  implementation  of  Model  fa  was  straightforward  but  nont rival  and  closely 
follows  that  of  Model  12.  The  implementation  of  Model  fa  was  limited  to  the  2D  and  3D 
solid  elements  and  it  is  compatible  with  all  of  the  other  options  in  ADIN  '.  However,  the 
user  is  again  warned  that  the  present  effort  validated  Model  fa  only  for  the  problems 
described  herein. 

Logic  was  modified  in  the  MAMLI2  (2D)  and  MATWUT  (3D)  subroutines  to  read 
and  print  the  new  material  property  information.  In  the  2D  and  3D  solid  element 
overlays  there  are  subroutines  present  for  each  material  model  (1-13).  Dummy 
subroutines  LLT2U6  and  F.  L I  3  Do  for  Model  fa  existed  and  these  were  replaced  with 
functional  logic  that  called  the  following  suite  of  new  subroutines  that  were  added; 


2D 

IT  H20fa 
T  1 1200 
M  R2  0b 
G  M  2  DO  fa 


31) 

II  If  306 
T  1 130b 
M  3  DOG 
G  M3  L)0fa 


While  the  corrections  were  straightforward  they  were  nontrivial  and  fairly  extensive  in 
that  approximately  1000  lines  of  code  were  added  or  modified. 

3.3  Data  Cards  for  MODEL  6 

Model  6  is  very  similar  to  Models  12  and  3  in  the  method  used  to  store  and 
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retrieve  material  properties  from  the  array  PUOP  as  a  function  of  strain  and 
temperature.  While  Model  b  can  he  used  without  strain  dependency  to  model 
temperature  dependent  anisotropic  materials,  this  usage  is  inefficient  and  not 
recommended.  Model  12  is  intended  for  these  problems  and  is  easier  to  use.  The 
material  properties  are  stored  in  a  table  in  a  manner  similar  to  the  pseudo  four¬ 
dimensional  array  I*  KOI*  described  on  the  succeeding  pages.  During  the  solution  phase 
the  properties  are  linearly  interpolated  as  illustrated  in  Figures  3-1  and  2-1.  Model  6  is 
biased  in  favor  of  strain  dependency  so  this  interpolation  is  performed  first.  If  more  than 
one  temperature  point  is  used  (M=NI*  AK(14)  >  1),  then  the  temperature  interpolation  is 
performed  on  the  strain  interpolated  properties.  Strain  and  temperature  points  must  be 

given  in  strict,  discrete,  ascending  order  ( e  <e  < .  .  .  <e ,  and  T  „  <  I  ,...<T  )  . 

12  k  12  m 

Note  that  the  properties  are  taken  as  undefined  outside  the  specified  strain  and 

temperature  range  and  execution  is  terminated  if  e  „  >  c  >  c ,  or  T  >  1  >  T 

1  k  1  m 

Note  also  that  the  strain  range  must  include  negative  values  for  compression  (see  Figure 
3-1).  Unlike  Model  12,  the  interpolation  for  Model  6  is  performed  on  the  data  for  the 
moduli,  Poisson's  ratio  and  thermal  expansion  coefficients.  These  interpolated  values  are 
then  used  to  define  the  strain-stress  coefficients  given  in  Figure  2-2  and  inverted  to  give 
the  stress-strain  law  given  in  Figure  2-3. 

The  strain  dependency  is  not  based  on  an  effective  strain,  because  as  described  in 
Appendix  C  this  would  be  a  nontrivial  computation  for  an  anisotropic  material,  hstead 
the  strain  dependency  is  of  the  following  form 


La 

= 

fca^> 

rb 

= 

l  b(Cb.  T) 

Ec 

= 

tc  T> 

Gab 

= 

Gab  ^ab'  7  ^ 

Gac 

= 

Gac  *Yac'  7) 

Gbc 

= 

Gbc  ^Ybc' 

vab 

= 

v><r> 

vac 

— 

vac  (T> 

vbc 

= 

vbc(r) 

ua 

= 

ua  (T) 

“b 

- 

%(n 

ac 

= 

nc(r) 

?A 


- - ^ 


Figure  3-1.  Model  6  Strain  Interpolation  Method. 
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I  his  approach  v/as  chosen  for  several  reasons.  I  he  moduli  are  usually  measured  in  a 
uniaxial  stress  test  and,  thus,  their  strain  dependency  is  primarily  related  to  the  uniaxial 
strain.  Data  for  the  strain  dependency  of  I'oisson's  ratio  is  rarely  available  and  when  it 
is,  it  is  usually  not  reliable.  Varying  the  I  I’oisson's  ratios  can  lead  to  numerical 
difficulties  that  may  not  be  easy  to  identify.  I  he  temperature  dependency  of  I’oisson's 
ratio  can  likewise  cause  problems  and  should  be  avoided  whenever  possible  (also  the  case 
for  Model  12).  However,  variation  of  the  uniaxial  and  shear  moduli  can  he  made  over  a 
large  range  without  numerical  problems.  The  coefficients  of  thermal  expansion  were 
taken  to  be  independent  of  strain  since  virtually  no  data  exists  to  identify  this 
dependency. 

A  small  number  of  modifications  to  the  ADINA  715  User's  Manual^  are  necessary 
in  order  to  describe  how  to  use  Model  6.  However,  there  are  no  changes  for  any  other 
material  models.  Kather  than  reproducing  extensive  portions  of  the  user's  manual,  in  the 
sequel  only  the  modifications  and  additions  to  the  manual  are  described. 

Modifications  for  the  2L>  and  3D  Solid  Element  Input  Data  (MODEL  6) 


Element  Group  Control  Card 
Columns  Variable  Entry 


45-48 


K=NI’  AR(I2)  Number  of  strain  points  ( » 2,  Model  6  only) 


53-56  M  — NP  AK(14)  Number  of  temperature  points  (Models  6  f*  12  only) 

57-60  NI’AR(15}~6  Material  model  number;  6  is  the  Orthotropic  Thermoelastic 

Strain  Dependent  Model  number. 


Set  NT  AK(1  5)--6  to  activate  tire  material  Model  6  input  data.  K“NP  AR(1 2)  defines 
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tf  number  of  str.iin  points  at  which  material  properties  ari1  to  bo  defined,  K  ^  2  and  is 
only  used  with  Model  (>.  M--NP  Ali  ('1 4)  defines  the  number  of  t  emperuture  points,  for  each 
strain  point,  for  whic  h  material  properties  are  to  he  input.  Nil*  AU(bl)  is  used  for  both 
Model  G  and  12,  but  with  Mode)  o  NPAK(/5J  can  be  given  as  1  (0  is  reset  to  1)  in  order 
that  strain  dependent,  temperature  independent  material  properties  can  be  defined. 


Columns 

Variable 

fcntry 

Card  Set  1-1 

1-10 

11-20 

• 

PKOP(1,N,1,1) 

P  K  00(2,  N, 1,1) 

• 

^all  rection  modulus  at  point  1,1 

C a2 1  a-direction  modulus  at  point  2,1 

• 

• 

Card  Set  1-2 

1-10 

11-20 

• 

• 

P  K  OP(  K,  N,1,1) 

1’ K  OP(  l,N,2,'l) 

P  R  0P(2,  N',2,1) 

• 

♦ 

• 

l^akl  a-£h rcct k>n  modulus  at  point  k,1 

1  a12  a-f<l rect'oa  modulus  at  point  1,2 

1  a22  a~c** r«ct ion  modulus  at  point  2,2 

• 

• 

• 

» 

P  l<  OP(  h,  N,2,1) 

• 

• 

t(1^2  a-direction  modulus  at  point  k,2 

• 

• 

Card  Set 1-M 

1-10 

11-20 

* 

PK0P(  l,N,M,1) 

P  K  0P(2,N,M  ,1) 

talm  a-direction  modulus  at  point  l,m 
t  a2m  a_d'r°Ction  modulus  at  point  2,m 

• 

• 

Card  Set  2-1 

1-10 

11-20 

• 

• 

• 

P  R  0P(  K,  N,  M  ,1) 

P  K  OP(  1,  N,1,2) 

P  K  OP£2,  N,  1,2) 

Eakm  a"^lr<JCt'on  modulus  at  point  k,m 

c"a11  a~direction  strain  at  point  1,1 

Ga21  a-d|rectlon  strain  at  point  2,1 

« 

• 

P  K  OP(  K,  N,  1,2) 

• 

Gak1  a_diroction  strain  at  point  k,1 
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Card  Set  2-2 


1-10 

11-20 

• 

P  KOP(  l,N,2,2) 
PROP(2,N,2,2) 

(  ;  a-diroLtion  strain  at  pf)irit  1,2 

>-a-,2  a-dirertion  strain  at  point  2,2 

• 

• 

P  R  OP(K,N,2,2) 

• 

• 

e  a-direction  Strain  at  point  k,2 

Card  Set  2-M 

1-10 

11-20 

% 

P  KOP(1,N,M,2) 

P  KOP(2,N,M,2) 

• 

t  ,  a-di rection  strain  at  point  1,m 
aim  *  ' 

e  2m  a-direction  strain  at  point  2,m 

• 

• 

• 

Card  Set  3-1 

1-10 

11-20 

• 

• 

PROP(K,N,M,2) 

P  K  OP(1,  N,1 , 3) 

P  R  OP(2,  N,1,3) 

• 

• 

• 

e  .  |n  a-direction  strain  at  point  k,m 

b-direction  modulus  at  point  1,1 
b-direction  modulus  at  point  2,1 

• 

• 

• 

• 

• 

PROP(K,N,1,3) 

• 

• 

^bkl  ^-direction  modulus  at  point  k,1 

• 

Card  Set  3-M 

1-10 

11-20 

• 

PKOP(1,N,M,3) 

PR  OP(2,  N,M,3) 

• 

b-direction  modulus  at  point  1,m 

d  1  m 

E^2m  b-direction  modulus  at  point  2,m 

• 

• 

• 

Card  Set  4-1 

1-10 

11-20 

• 

• 

P  R  O  P  (  K ,  N,  M  ,  3 ) 

PROP(1,N,1,4) 

P  KOP(2,N,1,4) 

• 

• 

• 

E b k m  b“c*'rcct'on  modulus  at  point  k,m 

b-direction  strain  at  point  1,1 
eb21  b-direction  strain  at  point  2,1 

• 

• 

• 

• 

• 

P  K  OP(  K,  N,1 ,4) 

• 

* 

ebk1  b-direction  strain  at  point  k,1 

Card  Set  4-M 

1-10  P R OP(1 , N,M ,4)  f-bim  b-direetion  strain  at  point  1,m 

11-20  P  K  OP(2,  N, M  ,4)  f  b-direction  strain  at  point  2,m 


,:bkin  b-direction  strain  at  point  k,m 
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P  K  OP(  K,  N,  M  ,4) 


Card  Set  5-1 


I- 10  I*  K  1)1*1 1,  N,  1,5) 

II- 21/  PKOP(2,N,1,r>) 


L  i  |  c-direction  modulus  at  point  1,1 
I  2i  c-direction  modulus  at  point  2,1 


I* K OP(  K,  N,  1,5)  I  t|_i  c-direction  modulus  at  point  k,1 


Card  Set  5-M 


1-10 

11-20 


I*  l<  OI’(  I ,  N,  M  ,5) 
F*  R  0IJ(2,  N,M  ,5) 


I  1(n  c-direction  modulus  at  point  1,m 
*'c2m  teetton  modulus  at  point  2,m 


Card  Set  6-1 


P  K  OI’(  K,  IM,  M  ,5)  *  ckm  c-(*' rect  ion  modulus  at  point  k,m 


I- 10  PROP(1,N,1,6) 

II- 20  P  K  0P(2,  N',1 ,6) 


cc^  c-direction  strain  at  point  1,1 
eC2-\  c-direction  strain  at  point  2,1 


P  K0P(K,N,1,(>)  ^ck  i  c-^'rect'on  sit-in  at  point  k,1 


Card  Set  6-M 


1-10 

11-20 


P  KOP(1,N,M,b) 
P  K0P(2,N,M  ,6) 


et1rn  c-direction  strain  at  point  1,m 
ct2,n  c-direction  strain  at  point  2,m 


Card  Set  7-1 


PKOP(K,N, M,b)  e  ^  c-direction  strain  at  point  k,m 


I- 10  PKOP(1,N,1,7) 

II- 20  PKOP(2,N,  1,7) 


^abll  a^-shcar  modulus  at  point  1,1 
C.^i  ab-sliear  modulus  at  point  2,1 


PK()P(K,N,1,7)  ^abkl  a':)~s*’ear  modulus  at  point  k,1 


Card  Set  7-M 


1-10 

11-20 


P  K  OP(  1 ,  N,M  ,7) 
P  KOP(2,N,M,7) 


^‘ablm  a^,_s*lear  modulus  at  point  1,m 
^'al)'»m  a^-shear  modulus  at  point  2,m 


P  K  OP( K,  N,  M  ,7)  ^'abkm  ,’i)"sl"-<'r  modulus  at  point  k,m 


i 


Card  Set  »-! 


l-iu 

1  1-20 

PkOP(  l,N,1,8) 

1’  K  OP(2,  N,1 ,8) 

Y  ^  l  ah-sliear  strain  at  point  1,1 
y  ,1,71  all-shear  strain  at  point  2,1 

• 

('  R  Ol'f  U,  N,  1 ,8) 

Y  b|-i  ah-shear  strain  at  point  k,1 

Card  Set  8-M 

1-10 

11-10 

PKOP(1,N,M,8) 

PROP(2,N,M,8) 

Yab1m  a^-shear  strain  at  point  1,m 
Yab2m  ab-s hea r  strain  at  point  2, in 

• 

Card  Set  9-1* 

1-10 

11-20 

• 

P  K  OP(  K,  N,M  ,8) 

PROP(1,N,1,9) 

P  KOP(2,N,1,9) 

• 

• 

• 

Y  ibkm  ab~shear  strain  at  point  k,m 

Cacn  ac-shear  modulus  at  point  1,1 
^“ac'M  ac~s^ear  modulus  at  point  2,1 

• 

• 

• 

• 

P  R  OP(  K,  N,1,9) 

• 

• 

^ackl  ac_shear  modulus  at  point  k,1 

• 

• 

Card  Set  9-M* 

1-10 

11-20 

• 

P  k  OP(1,  N,  M  ,9) 

P  R  OP(2,  N,  M  ,9) 

• 

Gacim  ac-s*'ear  modulus  at  point  1,m 
Car2m  ac-shear  modulus  at  point  2, in 

• 

• 

• 

Card  Set  10-1* 

1-10 

11-20 

• 

• 

• 

P  R  OP(  K,  N,  M  ,9) 

P  R  OP(  1,  N, 1,10) 

PR  OP(2,  N,1 ,1  0) 

• 

• 

(T'ackm  ac-shear  modulus  at  point  k,m 

Yacii  ac-shear  strain  at  point  1,1 
YaC2i  ac-shear  strain  at  point  2,1 

• 

• 

• 

P  R  OP(  K,  N,1 ,1  0) 

• 

• 

Y'ackl  ac-shear  strain  at  point  k,1 

Card  Sets  9,  10,  11  and  12  are  omitted  for  21)  Solid  I.  lemon  ts 


Card  Set  10-M* 

1-10 

1’  k  ()1’(1 ,  N,  M  ,1 0) 

Y  at:-sliear  strain  at  point  1,ni 

11-20 

P  R  OP(2,  N,M  ,1 0) 

^  ic2m  ac-sl>uar  strain  at  point  2,m 

• 

• 

P  R  OP(  K,  N,  M  ,1 0) 

Y  ac-shear  strain  at  point  k,m 

Card  Set  11-1* 

1-10 

PI<OP(1,N,1,1  1) 

*■‘1x11  bc-shear  modulus  at  point  1,1 

11-20 

• 

P  U  OP(2,  N,1 ,1 1) 

G|  tic-shear  modulus  at  point  2,1 

• 

• 

• 

P  R  OP(  K,  N,1 ,1 1 ) 

C^cki  hc-shear  modulus  at  point  k,1 

• 

• 

Card  Set  11-M* 

1-10 

P  K  OP(  1,  N,M  ,1 1) 

G^cin,  bc-shear  modulus  at  point  1,m 

11-20 

• 

P  K  0P(2,  N,  M  ,1 1) 

GjJC2m  bc-shear  modulus  at  point  2,rn 

• 

• 

Pk  OP(  K,  N,  M  ,1 1) 

• 

• 

G|lckni  bc-shear  modulus  at  point  k,m 

Card  Set  12-1* 

1-10 

P  k  O P ( 1 ,  N,1 ,1 2) 

Y|)Cn  bc-shear  strain  at  point  1,1 

11-20 

• 

PkOP(2,N,1,12) 

• 

Y|)c2i  bc-shear  strain  at  point  2,1 

♦ 

• 

• 

• 

* 

• 

P  k  OP(  K,  N,1 ,1 2) 

Ybck1  bc-shear  Strain  at  point  k,1 

• 

• 

Card  Set  12-M* 

1-10 

P  k  C)P(1 ,  N,M  ,1  2) 

Ybclm  bc-shear  strain  at  point  1,m 

11-20 

• 

P  k  OP(2,  N,  M  ,1  2) 

• 

Y b t' 2 m  bc-shear  strain  at  point  2,m 

• 

• 

• 

• 

P  k  OP(  K,  N,  M  ,1  2) 

• 

Ybcl  m  bc-shoa r  strain  at  point  k,m 

Card  Set  13 

1-10 

P  k  OP(1,  N,1 ,1  3) 

T ^  temperature  at  point  1 

11-20 

• 

P  k OP(1, N,2,1  3) 

• 

T  •,  temperature  at  point  2 

• 

• 

• 

P  K  OP(1,  N,  M  ,1  3) 

f  temperature  at  point  m 

*  Car  cl  Sets  9,  10,  11  and  1 2  are  omitted  for  21)  Solid  Elements. 


t 
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Card  Set  14 


1-10 

PKOP(  1,N,  1,14) 

vab1  a*>" Poisson's  ratio  at  (joint  1 

1  1-20 

1'  K  ()P(  1,  N,2,11) 

v  .  ratio  at  point  2 

• 

1'  K  OI’(  1 .  N,  M ,1  4) 

v  |)m  ab-Poisson's  ratio  at  point  m 

Card  Set  15 

1-10 

P  K  OI’(1 ,  N,  1 ,15) 

v.,  ..  ac-Poisson's  ratio  at  point  1 

11-20 

I'K  OP(  1/1,2,15) 

v  2  ac-PoissonN  ratio  at  point  2 

• 

P  K  OP(1,  N, M  ,1  5) 

v  ac-PoissonS  ratio  at  point  m 

Card  Set  16 

1-10 

P  KOP(l,  N,1 ,10) 

vhc1  l,c"l’uisson's  ratio  at  [joint  1 

11-20 

I'K  OP(1,  N',2, 1  0) 

vb(  2  *)C"l,,,'sson's  ratio  at  point  2 

• 

P  K  OP(1 ,  N,  M  ,1  6) 

vbcm  t,c-P0'*S0,,'s  ratio  at  point  m 

Card  Set  17 

1-10 

PK  OP(1,N,1,17) 

ad1  a-direction  uxpan.  cocf.  at  [joint  1 

11-20 

• 

P  K  0P(  1,  N,2,1 7) 

* 

ai2  r<-*ct ion  expan.  coef.  at  [joint  2 

• 

• 

• 

• 

P  R  OP(  1,  N,  M  ,1  7) 

• 

• 

<»  a-direction  expan.  coef.  at  point  rn 

Card  Set  18 

1-10 

P  K  ( )  l’(  1 ,  M,1 ,1 0) 

b-direction  expan.  coef.  at  point  1 

11-20 

• 

PKOP{  1,N,2,1R) 

• 

ub2  rect«°ri  expan.  coef.  at  point  2 

• 

• 

• 

• 

P  R  0P(1,  N,M  ,18) 

• 

ci^  b-direction  expan.  coef.  at  point  m 

Card  Set  19 

1-10 

P  K  OP(1,  N,1, 1  9) 

<.t c i  c-direction  expan.  coef.  at  point  1 

11-20 

• 

P  R OP(1,  N,2,19) 

• 

<lc2  c i r o c t i o n  cxpan.  coef.  at  point  2 

• 

• 

• 

PK  0P(1,N,M  ,19) 

u  m  c-direction  expan.  coef.  at  point  rn 

Card  Set  20 

1-10 

P  K  OP(1 ,  N,1 ,20) 

1  K  l  1  reference  stress  f  ree  temperature 
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<.4  Sample  Problems  lor  Model  0 

I  lie  i;-.*i)nu,t  r ies  and  boundary  conditions  of  the  214  and  41)  Model  (>  sample 
problems  are  the  same  as  the  Model  I 2  sample  problems  shown  in  I  igures  2-4  and  2-5  and 
described  in  Section  2.4.  I  he  21)  Model  (>  problem  corresponds  to  t her  file  K2I4j  in 
Appendix  A  and  the  314  problem  to  file  K.’Ds.  I  he  AI4IMA  7 it  data  for  the  214  and  4t4 
Model  0  sample  problems  are  given  in  I  aides  4.1  and  4.2,  respectively,  both  the  strain 
and  temperature  variations  were  used  in  the  problems,  and  a  very  interesting  feature  of 
this  data  is  that  the  strain  range  is  a  function  of  temperature.  A  tot  a  I  of  4  strain  poinis 
v\i’  AK(l2))  and  2  temperature  points  (NI'AKfM))  were  used,  fable  4.5  summarizes  the 
data  used  for  the  214  414  problems. 


fable  4.4  Material  Properties 


1  emperature 

(a,b,c) 

i 

a 

0 

-1 

0 

0 

0 

0 

0 

1 

0 

100 

-1 

10x1u° 

100 

0 

10x10° 

100 

1 

10x10° 

In  functional  form 

the  moduli 

va ry  as  follow 

for  MODI  1  6  Satnple  Problems 


_L_U 

_±c 

Cab 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

20x1  ()° 

1x1  ()f 

10.46x10° 

20x10° 

ixio1 

0 

20x10° 

1x1 0f 

Lb  =  1  6.46  x  10°  psi  (1/100)  f(  h,  l) 
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uf  —  -.1  (>t> 


I  bus. 


f  I  -=  t  ^  ur/20  -  - ii. ,i  x  10 


-3 


•llld 


-.003)  =  .01  7 


and 


E,  =  1.509  x  10°  psi 


flic  problem  is  again  one  of  plane  strain  (c  =  0)  in  which  the  in-plane  stresses  are 


statically  determinant 


°r=0a  2  0 


%  =  °b 


—  =  13333  psi 

The  strain-stress  law  becomes 


ro 

i 

-.1/1.509 

-  .  1 

/  o  1 

J  c r' 

Vj 

- . 1/1 . 50V 

-  .  1 

1/1  .  509 
-.15 

-.15 

1/2 

and  solving  for  o  and  c(|  gives 


x  10 


-6 


t>x  =  -4000  psi 
£(,  =  a.24  x  10 


-3 


So  after  2  iterations  the  solution  has  converged  to  an  accuracy  of  approximately  1%  for 
the  hoop  strain. 
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Table  3.1.  2D  “odel  6  Sample  Problem  -  Tile  R2l)b  Appendix  A. 


2  2 

1  2 

0.  G.  10.  1G. 

2  2 


0.  1.  10.  i. 


X 

1 

1 

0 

1 

1 

1 

1 

20.25 

X 

4 

1 

0 

0 

1 

1 

1 

20. 25 

5. 

5 

X 

44 

1 

0 

0 

1 

1 

1 
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Table  3.1  (cont'd)  2D  Model  6  Sample  Problem  -  File  R2D5  Appendix  A. 
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Table  3.2.  3D  Model  6  Sample  Problem  -  File  R3D5  Appendix  A. 
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4.  VbC  I  OK  1/ AT  ION  AN1>  OR  I  IM  1/  A I  ION  Ob  ALIINA  ON  lilt  ASC 


1.1  General 

I  lie  ASC  is  a  unique  fourth  generation  digital  computer  that  is  capable  of 
executing  single  floating  point  operations  at  speeds  comparable  to  typical  third 
generation  mainframe  computers  {e.g.(  CUC  P(>UO).  (lie  mode  of  computation  of  single 
floating  point  operations,  typically  these  floating  point  operations  are  intermixed  with 
integer  or  logical  operations,  is  called  the  scalar  mode  on  the  ASC.  I  lie  ASC  is  also 
capable  of  performing  vectorized  floating  point  operations  up  to  20  times  faster  than 
scalar  operations.  I  he  term  vectorization  on  the  ASC!  refers  to  the  'pipeline*  execution 
of  the  same  sequence  of  floating  point  operations.  The  ASC  has  the  capability  of 
simultaneously  performing  multiple  pipeline  operations.  Typically,  the  time  required  to 
initiate  and/or  terminate  the  pipeline  process  is  approximately  the  same  as  the  time  to 
perform  a  single,  scalar  operation,  i.e.,  about  20  clock  cycles.  However,  once  the 
pipeline  is  initiated,  subsequent  floating  point  operations  require  a  single  cycle.  This 
very  powerful  feature  of  the  ASC  was  expected  to  substantially  improve  the  efficiency 
of  AiJlNA,  and  a  portion  of  the  present  effort  was  directed  at  vectorizing  ADINA. 

A  second  feature  of  the  effort  described  in  this  section  is  optimization.  On  the 
ASC  there  are  two  compilers;  a  "fast"  compiler  that  does  little  optimization  of  the 
assembler  instructions,  and  a  "slow"  compiler  that  seeks  to  develop  more  efficient 
assembler  code  from  the  FORTRAN  instructions.  These  two  compilers  are  apparently 
similar  to  the  ORT— 1  and  OR  1=2  options  on  the  standard  CDC  TI  N  compiler,  both  the 
fast  and  slow  compilers  have  4  levels  of  optimization  that  can  be  selected.  In  most 
results  reported  herein,  the  default  level  3  was  used. 

4.2  Vectorization  Considerations 

Vector  i  zable  f  OR  IRAN  coding  is  typically  in  a  tight  DO  loop  of  the  following 

form 

DO  10  I--  1,1000 
A(  I)  =  A(l)  +  I5(  f|*  C(  I) 

10  CON  I  INU  l 


f 
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Hie  major  computational  effort  in  the  above  loop  is  the  floating  point  multiplication  and 
addition.  Such  loops  are  commonly  used  to  assess  the  computational  efficiency  of  codes, 
and  the  effort  in  the  above  loop  is  rated  at  1000  flops  (float in"  point  operations).^ 
However,  even  a  cursory  inspection  of  the  assembler  code  generated  by  the  TOKTKAN 
compiler  will  reveal  that  a  significant  number  of  overhead,  nonvectori/able  operations 
such  as  memory  fetches,  stores,  integer  addition  and  testing  must  also  be  performed. 
However,  fourth  generation  computers  such  as  the  ASC,  Cl<  AY,  SI  AK,  etc.  are  capable 
of  completing  the  1000  passes  through  the  loop  substantially  faster  than  1000  arbitrary, 
scalar  floating  point  operations. 

I  well  written  finite  element  code  such  as  ADINA,  the  bulk  of  the  CO  ( Central 
Processing)  time  is  consumed  in  the  formation  of  the  element  stiffness  matrices  and  in 
the  solution  of  the  assembled  simultaneous  equations.  Some  processing  time  must  also  be 
spent  performing  Input/Output  (I/O)  tasks  and  other  non-CI’  operations  such  as  reading 
and  writing  to  low  speed  disc  storage.  During  the  formation  of  the  element  matrices,  a 
substantial  number  of  computations  similar  to  the  above  loop  are  performed.  However, 
the  range  of  these  loops  is  typically  small  (h  integration  points,  8  shape  functions,  8 
nodes,  etc.)  Thus,  while  many  repetitive,  similar  calculations  are  made,  their  ranges  are 
small.  Vectorizat  ion  is  most  efficient  when  many  passes  (100)  are  made  through  each 
loop.  Tor  this  reason,  the  most  important  subroutine  for  vectorization  is  CO  ISO 
the  column  solver  routine. 

The  word  length  on  the  ASC  is  32  binary  bits,  the  standard  length  for  IBM  look- 
alike  computers.  This  length  is  marginally  adequate  for  explicit,  scalar  floating  point 
operations.  However,  at  least  40  bits  are  required  to  retain  sufficient  accuracy  for  large 
implicit  calculations  such  as  those  performed  by  ADINA.  Tor  this  reason,  all  the  floating 
point  variables  in  ADINA  were  made  double  precision,  64  bits.  This  has  some  adverse 
effects  on  optimization  and  vectorization. 

4.3  benchmark  Results 

The  important  computations  in  a  typical  column^  or  f  rental^  type  equation 
solver  can  be  represented  by  the  following  nested  loops 
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DO  20  N— 1 ,1 000 
SUM  =0.0 
DO  10  T-  1,1 0000 
SUM--SUM  + A(l)»  A(  1+  N) 

10  CONTINUE 

TO  1  =  10  T+SUM/ 10000. 

20  CON!  INU  L 

A  loop  structure  of  tliis  type  has  been  benchmarked  on  many  different  computers,  and 
has  been  found  to  be  a  very  good  benchmark  for  purposes  of  comparing  the  computational 
efficiency  of  different  computers.  Since  there  are  10  million  passes  through  the 
innermost  loop,  the  computational  effort  is  rated  at  10  Mflops  (Million  floating  point 
operations). 

On  the  ASC  the  same  loops  were  run  using  3  compiler  options:  (1)  using  the  fast 
compiler  with  the  default  level  of  "optimization*  and  no  vectorization  (vectorization 
cannot  be  performed  with  this  compiler);  (2)  using  the  slow  compiler  with  the  default 
level  of  optimization  and  no  vectorization;  and  (3)  using  the  slow  compiler  with  the 
default  level  of  optimization  but  giving  the  compiler  space  to  perform  the 
vectorization.  Although  it  is  possible  to  modify  the  TOR  THAN  coding  to  make  direct 
calls  to  vector  processing  routines,  this  was  not  done  in  any  of  the  work  reported  herein. 

Table  4.1  lists  in  ascending  order  the  central  processing  (CP)  time  in  seconds 
required  to  execute  the  10  Mflop  benchmark  case.  Unless  specifically  noted,  the 
compilation  was  done  using  the  best  applicable  level  of  optimization  and  the  arithmetic 
was  done  using  at  least  60  binary  bits,  b  some  cases  the  comparable  single  precision 
time  is  also  given. 

The  entries  in  Table  4.1  should  be  very  carefully  studied  for  they  convey 
significant  information  about  very  important  parameters.  The  vectorized  Cray  time  is 
excellent,  reflecting  a  computational  speed  of  15  x  10^  floating  point  operations  per 
second.  However,  scalar  arithmetic  on  the  Cray  is  noticably  slower  than  the  COC  170 
(same  as  a  CDC  7t>00).  The  range  of  times  in  the  table  vary  by  a  factor  of  almost  100. 


1  able  4.1.  Cl'  Kun  Times  for  ID  Mflop  benchmark 


Cornput  or 

Cl'  l  ime  (Sec) 

Remarks 

Cray-1 

.69 

Vectorized 

ASC 

.62 

Vectorized,  Single  free 

ASC 

3.2 

Vectorized 

COC  176 

4.1 

Cray-1 

4.9 

UCS  CUC  175 

5.7 

SUN  COC  750 

6.6 

CUC  170 

7.2 

No  Opt  (OPT— 1) 

liCS  COC  175 

7.3 

ASC 

16. 

Slow  Compiler,  Default  OPT 

C  DC  6600 

20. 

ASC 

26. 

Slow  Compiler,  No  OPT 

UN  IV  AC  iioo/ao 

34. 

Single  Prec 

UN  IV  AC  1100/00 

48. 

ASC 

56. 

Fast  Compiler,  Default  OPT 

The  ASC  is  seen  to  be  very  fast  for  single  precision  arithmetic,  but  double  precision 
increases  the  run  time  by  a  factor  of  4.  This  is  most  unfortunate  for  AL)  IN  A  (or  any  implicit 
code)  since  it  cannot  be  run  in  single  precision  with  accuracy.  However,  explicit  codes  such  as 
finite  difference  or  hydro-codes  that  can  use  single  precision  even  for  structural  response 
calculations  might  run  very  efficiently  on  the  ASC. 

The  importance  of  vectorizat ion  is  also  shown  in  the  table.  The  slow  compiler  ran  5 
times  faster  when  vectori/.ation  was  used.  This  is  a  big  improvement  but  substantially  less  than 
the  theoretical  maximum  factor  of  20.  Since  the  limit  on  the  inner  loop  was  10,000,  the 
overhead  associated  with  pipeline  initiation  should  be  negligible.  Thus,  it  would  appear  that  the 
best  vectorization  the  compiler  can  effect  for  double  precision  is  a  5  fold  reduction  in  run 
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t  i'lie.  Also,  tlu;  1  6  second  run  time  with  t  Ik;  slow  (opt  im  i  z  jug)  compiler  is  barely  faster  than  the 
CDC  6600.  I  his  implies  that  the  scalar  speed  of  the  ASC  in  double  precision  is  quite  slow.  The 
run  time  from  the  fast  compiler  seems  incredibly  slow,  even  slower  than  a  UN  IV  AC.  It  is  hard 
to  understand  a  factor  of  3.5  between  the  fast  and  slow  computers  on  the  ASC.  On  the  CUC 
1'/t>  the  factor  is  less  than  2. 

AO  IN  A  benchmark 

While  the  loop  structure  accurately  reflects  the  equation  solver  effort,  the  final 
f  evaluation  must  be  of  actual  runs  o..  AOINA.  I  he  table  below  lists  the  ASC  and  CDC  17b 

;  CT  times  for  six  20  and  six  30  runs.  All  of  the  ASC  runs  used  the  best  optimized  and 

1  vectori zed  version  of  AOINA  that  was  developed  during  this  effort. 

I 

j  Table  4.2  ASC  vs.  COC  AOINA  Run  Times  (CT*  seconds) 


2  0  3  0 


ASC 

C  PCI  76 

Ratio 

_ASC_ 

COC176 

Ratio 

1  .81 

.22 

3.68 

9.0 

1.90 

4.73 

2  1.01 

.23 

4.39 

8.86 

2.01 

4.40 

3  1.14 

• 

ce 

4.07 

1.21 

.37 

3.27 

4  4.58 

.73 

6.27 

18.56 

6.79 

2.7  3 

5  1.80 

.39 

4.61 

10.62 

2.25 

4.72 

6  13.97 

1.89 

7.59 

53.81 

9.16 

5.87 

Average 

5.1 

Average 

4.3 

The  ASC  is  seen  to  run  AOINA  slower  than  the  COC  17b,  typically  4-5  times 
slower,  h  one  30  case  the  factor  was  less  than  3,  and  the  efficiency  on  the  ASC  seemed 
to  improve  as  the  problem  size  increased.  The  above  problems  are  obviously  small  with 
run  times  not  exceeding  1  minute  on  the  ASC.  The  small  size  might  bias  the  results 
unfavorably  against  the  ASC. 

Referring  to  Table  4.1,  the  COC  17b  was  found  to  be  approximately  4  times  faster 
than  the  optimized  ASC  running  in  a  scalar  mode  (no  vectoriz.ation).  1  he  ratios  in  lable 
4.2  suggest  that  the  vector  ization  of  AOINA  may  have  had  little  effect  on  the  run  time. 
This  suggestion  is  supported  by  the  run  times  recorded  compiling  AOINA  with  the  fast 
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compiler  (these  results  have  not  been  reported)  which  generally  run  only  2  times  slowei 
than  wit  It  the  best  opt  i  mi  /  at  ion  and  vec  tori  /at  ion.  Agai  n  referring  back  to  I  able  4.1,  the 
fast  compiler  ran  3.5  times  slower  than  the  scalar  mode  of  the  optimizing  compiler  and 
17.5  times  slower  than  the  best  optimized  and  vectorized  compiler.  However,  the  10 
Mtlop  benchmark  is  strongly  dominated  by  pipeline  floating  point  operations,  and  in 
ADINA  nearly  W of  the  executable  statements  involve  non-floating  point  or  scalar 
floating  point  operations.  finis,  it  is  not  unreasonable  that  the  run  times  with  ADIN\ 
would  be  slower  in  relation  to  the  OK  17b  (a  very  fast  scalar  computer)  than  the  10 
M  flop  benchmark  or  that  the  slow  compiler  would  be  less  effect  ive  on  AD  IN  A  than  on  the 
benchmark* 

4.4  Discussion 

The  results  of  the  vectorizat ion  and  optimization  studies  presented  in  this  section 
are  surprising  and  disappointing.  I  he  best  ASH  double  precison  run  of  the  10  Mflop 
benchmark  was  20&  faster  than  a  CDC  176,  but  the  best  ASC  run  times  on  ADINA  are  4 
times  slower  than  the  CDC  17b.  Without  further,  more  detailed  studies  it  is  not  possible 
to  be  certain  about  the  cause  of  the  ADINA  inefficiency.  However,  the  most  likely  cause 
is  the  very  high  percentage  of  scalar  operations  in  ADINA  and  the  very  poor  efficiency  of 
the  AS C  in  performing  scalar  computations.  It  is  unlikely  that  further  significant 
improvements  can  be  made  in  ADINA  using  the  current  best  (slow)  compiler.  It'  many 
large  problems  (30  minutes)  are  to  be  run,  it  would  be  worthwhile  to  force  vectorization 
in  the  equation  solver  by  modifying  the  FORTRAN.  The  bulk  of  ADIN  A  should  be  lef  t  as 
is  and  future  improvements  in  efficiency  should  he  done  using  small  subsets  of  the 
ADINA  code,  e.g.,  equation  solver,  assembler,  element  stiffness  routines,  etc. 

At  the  start  of  this  effort,  we  were  told  that  the  slow  compiler  was  unreliable  and 
often  ran  5  times  slower  than  the  fast  compiler.  Our  findings  were  that  the  slow 
compiler  was  quite  reliable  and  that  it  typically  ran  3  times  slower.  The  CDC  FIN 
0PT=2  compiler  typically  runs  2  times  slower  than  Of’T=1  and  its  object  code  typically 
runs  2  times  faster.  We  found  no  reason  to  use  the  fast  compiler  and  recommend  that  its 
use  be  discontinued. 

The  main  reason  that  the  ADINA  run  times  are  so  slow  is  probably  due  to  the 
compiler.  The  run  times  on  the  10  Mflop  benchmark  varied  from  .fi  to  50  seconds.  This 
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is  an  unacceptable  range.  Reliable  compiler*  that  never  take  longer  than  10  seconds  in 
single  precision  should  be  provided.  Also,  it  is  surprising  that  scalar  operations  are  so 
slow.  Ag.ain,  we  suspect  the  compiler  or  operating  system.  Our  opinion  is  that  until 
basic  improvements  in  the  ASC  are  made,  ADINA  will  continue  to  run  4  times  slower 
t  han  on  a  C  l)C  1  7b. 
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5.  CONCLUSIONS 


The  work  reported  herein  lias  lead  to  significant  improvements  in  the  AO  IN  A  code* 
on  NR  L's  ASC  computer.  Material  Model  12  was  added  to  the  20  and  '3  0  solid  element 
types  in  order  to  permit  the  modeling  of  anisotropic  thermoclustic  materials,  and 
mate  ial  Model  6  was  added  to  these  same  element  types  to  permit  the  modeling  of 
strain  dependent  nonlinear  anisotropic  thermoelastic  materials.  These  additions  should 
markedly  improve  NR  L's  ability  to  accurately  determine  the  response  of  naval  structures 
to  intense  rapid  heating. 

Considerable  effort  was  also  directed  at  improving  the  efficiency  of  AUINA  on 
the  ASC.  Through  use  of  optimization  and  vectorization,  AUINA  now  runs  at  least  twice 
as  fast  on  the  ASC.  Unfortunately,  AUINA  still  runs  four  times  slower  on  the  ASC  than 
on  a  CUC  17b.  As  discussed  in  Section  4,  this  is  probably  caused  by  the  poor  compiler 
provided  by  FL  I  he  current  work  did  not  attempt  the  solution  of  any  large  scale 
problems  and  it  is  likely  that  the  ASC  will  improve  its  ratio  to  the  CUC  176  on  large 
problems.  It  is  possible  that  by  modifying  the  FOR T  AN  instructions  in  ADINA  to  make 
direct  calls  to  assembler  language  vector  processing  routines,  the  run  times  can  be 
further  reduced.  However,  it  is  doubtful  the  AUINA  will  ever  run  faster  on  the  ASC  than 
on  a  COC  176. 

A  great  deal  has  been  done  to  improve  ADINA,  but  more  remains  to  be  done.  1  he 
next  Section  6  gives  a  detailed  list  of  recommendations  for  further  work  in  this  area. 


6.  RECOMMENDATIONS 


Much  has  been  dune*  to  make  AD  IN  A  operational  and  run  efficiently  on  the  ASC, 
but  much  remains  to  be  done.  The  following  is  a  list  of  recoin  mended  improvements  and 
modifications  to  AD  IN  A. 

6.1  Development  of  a  Preprocessor 

ADINA  has  no  preprocessing  or  data  generation  capabilities  and  its  input  structure 
is  exceeding  complex,  occassional ly  redundant,  and  inconvenient  to  use.  This  condition 
detracts  from  effective  utilization  and  wastes  many  hours  of  labor  in  the  data 
preparation  process.  A  preprocessor  with  some  form  of  data  generation  should  be 
developed  as  soon  as  possible.  At  a  minimum  this  should  include  a  free  format  reader 
and  nodal  point  and  element  generator  such  as  T  t  XG  AP-2  0^ ^  or  1  E  XC  AP-3  0^ 

6.2  Development  of  a  Postprocessor 

The  output  from  ADINA  consists  entirely  of  printed  output,  no  graphical  output  is 
available.  However,  ADINA  does  have  ‘portholes'  that  permit  some  user  interaction 
with  the  code  (for  both  pre  and  post  processing).  A  postprocessing  module  should  be 
developed  to  identify  maximum  and  minimum  stresses  and  strains  by  material  and 
element  group  and  to  plot  element  connectivities,  deformed  grid,  and  stress  and  strain 
contours  on  parametric  surfaces.  This  should  be  done  using  a  ‘neutral  plot  file"  concept 
such  as  used  by  Pacifica  Technology  intheTEXGAP  codes. 

6.3  Render  all  of  ADINA  70  Operational  on  the  ASC 

At  present  only  the  20  and  3  0  solid  element  types  are  operational  for  static  and 
transient  analyses.  The  other  element  types  such  as  truss,  beam,  shell  and  fluid  elements 
should  be  made  operational  and  fully  checked  out.  Several  auxiliary  options  such  as 
frequencies,  creep,  etc.  also  should  be  made  operational. 
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t>.4  Further  Efficiency  Improvement 


A  low  level  of  effort  as  described  in  Section  4.4  should  be  continued  in  Older  to 
.jsc  .-rtain  the  most  efficient  method  of  running  ADIS' Aon  the  ASC. 

b.'->  Local  Orthotropic  Material  Properties 

Most  of  t  fie  anisotropic  materials  that  are  used  in  typical  structures  are 
transversely  or  locally  orthotropic.  The  local  orthotropic  axes  usually  vary  as,  for 
example,  with  axisymmetric  shells,  fhe  capability  to  model  these  locally  orthotropic 
materials  and  variable  local  axes  orientation  should  be  added  to  ADINA  to  permit  the 
more  effective  modeling  of  these  materials. 

6.6  Development  of  a  Mew  Shell  Element 

The  shell  element  implemented  in  ADINA  is  suitable  mainly  for  moderately  thick 
shells  and  transitions  from  shell  to  solid  (continuum)  element  types.  This  element  is 
unsuitable  for  general  shells  and  it  is  particularly  inefficient  for  thin  shells,  a  class  of 
problems  of  considerable  interest  to  t fie  Navy.  A  new  shell  element  should  be  developed 
and  added  to  ADINA,  retaining  the  current  element  as  an  option.  A  stiffener  (beam) 
clement  that  is  compatible  with  the  shell  element  should  also  be  developed  since  many 
naval  shell  structures  are  stiffened. 

t>.7  Addition  of  Material  Models  to  the  Shell  Element 

Presently  the  shell  element  is  limited  to  isotropic  elastic  or  plastic  material 
models.  No  thermoelastic  or  ortliotropic  materials  can  be  modeled,  fhe  MODE  l.=2,3,6  & 
12  material  models  should  be  added  to  the  shell  element. 


6.8  ADINA-T  Interface 

NK  L  has  put  considerable  effort  into  implementing  and  optimizing  both  ADINA 
and  ADINA-T  on  the  ASC.  The  interface  between  these  codes  should  be  thoroughly 
tested  and  upgraded  to  meet  the  needs  of  NK  I.. 
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O.'t  fc valuation  &  Improvement  of  Nonlinear  Algorithms 


As  described  in  Appendix  II,  there  are  open  questions  about  the  nonlinear  solution 
algorithms  used  in  ADINA.  1  he  code  has  been  thoroughly  checked  out  on  linear  static 
and  dynamic  and  static  nonlinear  elastic  problems.  However,  even  highly  skilled  users 
have  reported  problems  with  AD  IN  A  most  notably  for  dynamic  plasticity  problems.  The 
present  version  of  ADINA  does  not  permit  variable  (adaptive)  time  steps.  Ibis  is  a 
serious  deficiency  that  must  be  removed  if  dynamic  nonlinear  problems  are  to  be  run 
efficient  ly. 

6.1 0  Improved  Documentation 

While  the  documentation  of  ADINA  is  voluminous,  it  is  severely  limited  in  scope, 
fails  to  give  the  user  necessary  guidance  or  recommendations,  and  is  very  difficult  for 
thi;  uninitiated  to  read  and  learn  from.  Personnel  at  the  NR  l.  and  other  Navy  labs  would 
be  able*  to  make  much  more  effective  and  efficient  use  of  the  code  if  better 
documentation  were  available. 

ADINA  is  coded  in  easy-to-read  FORTRAN,  and  it  is  organized  into  modular, 
functional  overlays  that  perform  identifiable,  logical  tasks.  Thus,  ADINA  is  relatively 
easy  to  modify.  For  example,  .Hiding  a  new  element  type  such  as  an  improved  shell 
element  or  adding  new  material  models  such  as  those  described  in  Sections  2  and  3  can  be 
done  by  a  skilled  programmer  in  a  reasonable  period  of  time.  The  personnel  using  ADINA 
could  readily  learn  bow  to  make  these  modifications  themselves  if  a  Programmer's 
Manual  was  available.  A  Programmer's  Manual  typically  describes  the  organization  of 
the  code  and  the  tasks  performed  by  individual  routines  or  groups  of  routines  as  well  as 
describing  the  COMMON  blocks,  addressing  and  Input/Output  structure.  We  recommend 
that  such  a  document  be  prepared. 
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APPENDIX  A 

ADINA  FILES  ON  IMF  ASC 


This  appendix  describes  t lie 
use.  Any  questions  about  these 
I  echnology ,  11696  Sorrento  Valley 


ADINA  files  that  have  been  left  on  the  ASC  for  future 
files  should  be  directed  to  Dale  Ranta  at  Pacifica 
Road,  San  Diego,  California  92121,  (714)  453-2530. 
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U  StUCA I  / 1  J(.  y  I  i»0/  IJU  N  H  IM 


Most  of  the  ADINA  files  are  on  this  path 
M ACKOS 

M  AC K O  to  Get  2-0  ADINA  files  from  tape 
M  ACKU  to  Get  3-D  AD  IN  A  files  from  tape 
M  ACKO  to  Get  2-O/S-O  ADINA  files  from  tape 
10  Mflops  lienchmark  problem 


Patl)  Name 
Sons: 

NAM  LS 
A  M  A  C 
L  IliU 
LIB 

AL  INK 
ABS2L) 
XQT2D 
ABS3D 
X  QT  3  O 
L  INK2D 
L  INK3D 
SOUK  Cf 
UP 

U  P  L)  A  T  E 
RUNS 
OU1 PUT 


USE  K  C  AI/D(.3/Bau/DUNHK  1/AUtNA 

Names  of  ADINA  files 
M  ACKOS 

C  II  E  K  source  I  ibrarys  on  t  his  path 
C  IF  E  K  execution  library  s  on  this  path 
M  AK  CO  to  link  2-D/3-U  ADINA 
Absolute  element  for  2-1)  ADINA 
M  AK  CO  to  execute  2-i)  AUINA 
Absolute  element  for  3-0  AUINA 
M  AK  CO  to  execute  3-D  AUINA 
M  AK  CO  to  create  ABS2D 
MACRO  to  create  ABS3U 
Source  code  for  M0  0LL12 
M  Alt  COS 

C  IF  E  R  update  decks  are  on  this  path 
ADINA  test  problems  are  on  this  path 
ADIN  A  test  problem  output  is  on  this  path 
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'I • '} N‘1,1U‘  USE  K  C  A3 /Dbyttao/DUNII K  v  l  iuu 

SONS: 


AO  IN  A 

- 

C  II  1:  K 

source 

for  the  ‘main*  AO  IN  A  overlay 

AO  IN  1 

- 

Cll  1  11 

source 

for  the  "input’  ADINAoverlay 

1  DOM  1  1: 

- 

C  II  1. 11 

source 

for  the  "input"  ADINAoverlay 

l.  0 1  2  1  J.i 

- 

Clb  I  II 

source 

for  the  'input’  ADINAoverlay 

L  0 1  2  04 

- 

ClFbll 

source 

for  the  "input"  AOINA  overlay 

1  Ol  200 

- 

Cll  b  11 

source 

for  the  "input"  ADINAoverlay 

b  Ol  2  07 

- 

C  IF  L  11 

source 

for  the  "input"  AOINA  overlay 

i  oi  20a 

- 

C  IF  b  11 

source 

for  the  "input*  ADINAoverlay 

t  02  OIU 

- 

C  IFF  11 

source 

for  the  "input"  AOINA  overlay 

1:  02  O  l  2 

- 

CIIER 

source 

for  the  "input"  ADINAoverlay 

b  02  013 

- 

CIFEK 

source 

for  the  "input"  AOINA  overlay 

F  02  014 

- 

Cll  Ell 

source 

for  the  “input"  ADINAoverlay 

1  1 1 II  b  DM 

- 

ClFbll 

source 

for  the  "input"  AOINA  overlay 

1  013  03 

- 

CIFEK 

source 

for  the  "input"  ADINAoverlay 

L  OF 3  04 

- 

ClFbll 

source 

for  the  "input"  AOINA  overlay 

1:  OT  3  06 

- 

CIFEK 

source 

for  the  "input"  ADINAoverlay 

1:  o  r  3  07 

- 

C  IFE  K 

source 

for  the  "input"  Ab)  IN  A  overlay 

E  LJ 1  3  DU 

- 

CIFEK 

source 

for  the  "input"  ADINAoverlay 

I.  03  010 

- 

ClFbll 

source 

for  the  "input"  AOINA  overlay 

b  03  012 

- 

ClFbll 

source 

for  the  "input*  ADINAoverlay 

1  /U  M  M  Y 

- 

ClFbll 

source 

for  the  "input"  ADINAoverlay 

13  LOCK 

- 

ClFbll 

source 

for  the  "input"  ADINAoverlay 

1  O  A  0 

- 

CIFEK 

source 

for  the  "load"  ADINAoverlay 

Path  Name 


US  I:  K  C  AI/D03/li80/DUNIIK1/b  111 


SONS: 

AD  IN  A 

- 

C  IFER 

execution 

1  ibrary 

tor 

t  be 

'main'  AU  IN  A  overlay 

AO  INI 

- 

C  lb  I  K 

execution 

library 

for 

the 

’input’  ADINAoverlay 

1  ODM  EE 

- 

CIFEK 

execution 

library 

for 

t  he 

'input'  ADIN  A  overlay 

I:  o  r  2  03 

- 

C  II  I  K 

execution 

library 

for 

the 

’input’  ADINAoverlay 

E  O  l  2  04 

- 

CIFEK 

execution 

library 

for 

the 

‘input’  ADIN  A  overlay 

E  Ol  206 

- 

C  IF  EK 

execution 

library 

for 

the 

’input’  ADINAoverlay 

t  Ol  2  07 

- 

CllbK 

execution 

1  ibrary 

for 

t  he 

’input’  ADINAoverlay 

t  O  l  2  08 

- 

CIFEK 

execution 

library 

for 

the 

’input’  ADINAoverlay 

1  02  DIO 

- 

C  lb  b  K 

execut  ion 

1 ibrary 

for 

the 

’input’  ADINAoverlay 

t  02  012 

- 

ClFbK 

execution 

library 

for 

the 

’input*  ADINAoverlay 

l  02  013 

- 

C  IF  b  K 

execut  ion 

1 ibrary 

for 

t  he 

’input’  ADINAoverlay 

E  02  014 

- 

CIEEK 

execution 

library 

for 

the 

'input*  ADINAoverlay 

T  I  I  K  b  OM 

- 

ClFbK 

execution 

1 ibrary 

for 

the 

’input*  ADINAoverlay 

b  01  3  03 

- 

CllbK 

execution 

library 

for 

the 

’input’  ADINAoverlay 

E  O  l  3  04 

- 

CIFEK 

execution 

library 

for 

t  he 

’input’  ADINAoverlay 

E  O  r  3  06 

- 

CIFEK 

execution 

library 

for 

the 

’input’  ADINAoverlay 

E  OT  3  07 

- 

CIFEK 

execution 

library 

for 

the 

‘input’  ADINAoverlay 

E  OT  3  08 

- 

CIFEK 

execution 

library 

for 

the 

’input"  ADINAoverlay 

F.  03  010 

- 

CIFEK 

execution 

library 

for 

the 

’input’  ADINAoverlay 

E  03  01  2 

- 

CIFEK 

execution 

library 

for 

the 

‘input’  ADINAoverlay 

OUM  M  Y 

- 

CIFEK 

execution 

1 ibrary 

for 

the 

’input’  ADINAoverlay 

13  I.OCK 

- 

CIFEK 

execution 

library 

for 

the 

’input’  ADINAoverlay 

bO  AO 

- 

C  IF  E  K 

execution 

1 ibrary 

for 

the 

’load’  ADINAoverlay 

Pat  h  Name 


UStKCAT  /  Ub  3/13  «0/  UU  (si  1 1 K1/KU  N  S 


SUNS: 

K  2  Ul 

- 

AUINA  2-U  Test  problem  1 

u  2 1  n 

- 

AUIN  A  2-U  1  est  problem  2 

K2  U3* 

- 

AUINA  2-~  lest  problem  3 

K  2  U4 

- 

ADINA2-U  Vest  problem  4 

R 2  US* 

- 

AUINA  2-L>  lest  problem  S 

K  2  i  )b 

- 

Al) IN  A  2-U  1  est  problem  b 

K  3U1 

- 

AUINA  3-U  Test  problem  1 

K  3  U2* 

- 

AUIN  A  3- U  l  est  problem  2 

R3U3 

- 

AUINA  3-U  Test  problem  3 

l\  3  134 

- 

AUINA  3-U  Test  problem  4 

R  3  US* 

- 

AUINA  3-U  Test  problem  5 

K  3  Ut> 

- 

AUINA  3-U  Test  problem  b 

•Sample  problems  described 


in  Sections 


and  3 . 
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APl'tNDIX  (5 
fcXf  CUI  ION  I’KOULI  M  S 


I'criiintK1 1  lit  Pacifica  I  ei  hnnlogy  have  been  using  the  AO|\A  rode  since  its  first 
release  in  10/ j.  While  we  have  found  it  to  he  a  vert  reliable  code  that  solves  a  very 
large  class  oi  problems,  we  have  also  observed  sarious  execution  errors.  Almost  all  the 
difficulties  we  have  noted  are  associated  with  tee  plasticity  models  and  to  a  lesser 
degree,  the  constant  time  step.*  this  is  especially  true  for  plasticity  problems  with  very 
rapid  loading.  AINNA  almost  always  fails  on  tins  t',pe  of  problem,  except  when  the  time 
step  si/e  is  severely  reduced  (usually  less  than  the  Courant  stability  limit  At  c  Ax/c  ). 

Figures  b-1  and  t j- 2  illustrate  a  typical  dynamic  plasticity  problem  that  ADlNAhas 
been  unable  to  solve.  An  excellent  numerical  solution  to  this  problem  was  computed  with 
the  UY  NAP  LAY  code,^J^  and  recently  a  solution  was  also  achieved  with  the  Si  ACS-C1 
code  (this  geometry  was  also  used  in  an  experiment  reported  by  Wu  and  Witmcr).  Our 
AUINA  run  used  a  5,, sec  time  step  (certainly  too  large'  for  good  accuracy)  and  at  the  end 
of  the  first  time  step  in  the  element  at  the  cm.ui  the  stress  was  ~07,000  psi  and  the 
radial  displacement  was~.023".  I  bus,  the  strain  from  the  strain  displacement  equations 
was  l  -  u/r  -.OOtf,  but  from  the  stress-strain  law  the  plastic  strain  was  c  "((i7,(.KK>- 
42,(IUU)/.7ti7  x  103-  .32,  about  40  times  the  strain  consistent  with  the  displ acement s.  In 
our  opinion,  this  error  is  caused  by  a  lack  of  closure  in  the  plasticity  routines,  that  is,  a 
failure  to  test  the  stress  state  again* t  the  c.omput ed  strain  state. 

These  difficulties  are  potentially  serios.w  and  warrant  further  invest  igat  ion. 
However,  NR  L  is  not  presently  using  AUINA  for  this  type  of  problem.  It  is  possible  that 
this  deficiency  can  be  easily  fixed,  but  this  will  not  be  known  until  a  cure  is  attempted. 


*  It  is  remarkable  that  bathe  has  been  able  to  retain  the  constant 
time  step  in  a  general  purpose  nonlinear  code. 
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Figure  B-l.  Typical  Dynamic  Plasticity  Problem  for  ADINA. 


Figure  B-2.  Stress-Strain  Law  for  Plasticity  Problem. 


APIM  .NOIX  C 

AN  ISO  r  KOI’  1C  PI  FORM  Al  ION  IHLOKY  P  L  AS  I  IC  II  Y 


Refot mation  Iht-ory  I'lasticity  is  a  theory  that  approximates  the  nonlinear  stress- 
strain  relationships  of  plasticity,  hot  avoids  the  complexity  of  solving  the  actual 
plasticity  equations.  Within  reasonable  limits,  this  is  generally  a  good  approximation  of 
material  nonlinearity  for  loading  histories  without  unloading  and  where  the  ratio  of  the 
stress  components  do  not  significantly  change.  Reformation  theory  gives  results  that  arc- 
path  independent  and  the  theory  is  implemented  using  a  secant  modulus.  Deformation 
theory  can  be  readily  incorporated  into  the  ARINA  computer  code  because  although 
ARINA  performs  "incremental*  solutions,  it  uses  a  secant  modulus  stiffness  matrix  for 
the  Model  3  and  the  recently  developed  Model  12  temperature  dependent  properties,  and 
the  results  are  likewise  load  path  independent. 

The  following  development  is  for  anisotropic  deformation  theory  plasticity  and 
illustrates  how  it  could  he  incorporated  into  ARINA  in  a  manner  similar  to  Model  6. 
botropic  materials  would  follow  in  a  similar  manner,  and  note  that  while  the  material  is 
elastically  isotropic,  the  secant  (or  tangent)  stress-strain  relationship  becomes 
anisotropic  after  yield. 

Anisotropic  Stress-Strain  Law 

Derivations  are  limited  to  anisotropic  materials  that  have  stress-strain  laws  in  the 
elastic  range  that  are  expressed  in  the  principal  material  axes  as 

{of  =  [CCJ  (c-uATf  (C-1) 


where 


Hit  ii  res  Cla-CIf  i  I  lust  i  ate  the  uniaxial  (stress)  stress-strain  behavior  of  the 
material  in  the  piincipal  material  directions.  for  simplicity  attention  is  limited  to 
materials  that  have  a  linear  hardening  law.  lo  fully  describe  a  deformation  theory 
anisotropic  material  it  is  necessary  to  measure  '>  elastic  moduli  C.^,  C.|2»  V  ^21' 
l.'2H,  CiV  C44,  CV)(  Cb(>  (or  equivalently  the  l»  engineering  elastic  properties  tx,  t  .  L  ?, 
vx  ,  v  ,  v/x<  Cxy,  Gyz,  G/x),  c.  coelf icients  of  thermal  expansion  nx,  a  ,  u?,  n  ,  ayz, 
u/x,  b  yield  stresses  Yx,  Yy>  Y/(  Y  ,  Yy/,  Y/x,  and  b  hardening  moduli  1 1  x,  Hy,  H/f  Hxy( 
II  ,  HzX*  Thus,  there  are  27  material  properties  to  be  measured.  Also,  all  27  of  these 
material  properties  can  be  temperature  dependent  (the  9  elastic  moduli  and  6 
coefficients  of  thermal  expansion  are  functions  of  temperature  in  Model  12). 

Anisotropic  tffective  Stress 


following  closely  frank  Weiler's  development  of  the  anisotropic:  deformation  theory 
in  the  1)0  AS  IS  corle,^  the  existence  of  an  effective  stress  measure  (or  equivalently  an 
anisotropic  yield  function)  in  the  following  form  is  postulated 


a  =  ft 


) 


o  o  2 

1(  r  r 

x  y 

„  2 

ft  -i  +  ft  t  + 
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where  the  r  weighting  factors  are  chosen  lo  make  the  effective  stress  independent  of 
volumetric  strains,  thus,  making  the  effective  stress  a  measure  only  of  deformation 
producing  strain  states. 


rx  =  (C11  +  CVJ  *  C13)/C,  ry=(C12  ♦  C22  +  C23)/C,rz=(C13  +  *  CJ3)/C  (C-3) 

where  for  convenience  we  choose  C  =  +  C^2  +  C^3  <ind,  thus,  rx  -  1.  Now  define 

1!x  ••=  2i +  h 

fty  =  ft1+ft2  ( C-4a) 

ftz  -  ft2  -  ii3 


and  the  inverse  relations 
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Anisotropic  tffective  Stress- 

Strain  Relation 

Next  an  anisot  ropic  ef  feet  ive  st  ress-st  rain  rclat  ioriship  t  hat  is  the  same  form  as  the 
uniaxial  components  is  postulated,  This  is  illustrated  in  Figure  C2.  1  he  effective  yield 
stress  is  V  ,  The  elastic  modulus,  H  ,  I  he  hardening  modulus  is  I!  and 


1  =  1 

K  ft 

lor  convenience  the  effective  stress-strain  relation  is  chosen  to  he  the  same  as  the 
uniaxial  ov  -  r:  relation,  thus  Y  =  Y  ,  L  —  I'  ,  IT  -  H  and  e  =  t:  .  This  is 
arbitrary  because  the  effective  relation  must  embody  all  the  component  relations. 


(06) 


The  fi's  are  determined  by  equating  the  plastic  work  of  each  component  to  the 
plastic  work  of  the  ef  feet  ive  stress.  Without  reproducing  the  algebraic  detai  Is,  this  gives 
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Thus,  we  have  the  rather  unpleasant  result  that  the  , 's  depend  on  the  effective  stress  for 
an  anisotropic  material.  (The  's  are  constants  for  an  isotropic  material.)  Note  that 
if  -1  because  of  how  we  chose  the  normalizing  factor  for  the  r's  and  the  effective  stress- 
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st  rain  relation. 

Anisotropic  Deformation  I  lieory  flow  Kule 

Attain  omitting  some  algebraic  details,  the  total  plastic  strain  components  are 
expressed 


cf  -  A  S.  =  7-  S.  ,  i-x  ,  y  ,  7 . ,  xy  ,  y  z  ,  zx  (C-«) 

o 

where  A  is  a  positive  proportionality  factor,  t:  is  the  total  effective  plastic  strain  and 
-S-  are  the  so-called  anisotropic  deviatoric  stresses 
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fcffective  Secant  Anisotropic  Deformation  Theory  Stress-Strain  Relation 

The  secant  stress-strain  law  for  the  anisotropic  Information  theory  plasticity  is 
expressed 

t'fi  =ICC,,1  Is -a  A  1}  ( C-1  2a) 


whore 


lCep]  -  HL>eJ  +  A l Li  1 1 


( C-1 2h) 


anti 


lDej-lCe]"1  (C-12c) 

l fi ]  is  the  anisotropic  deviatoric  stress  matrix  given  in  (C-9)  and  [CR)  is  the  elastic  stress- 
strain  matrix  given  in  (C-1). 

Implementation  of  the  Anisotropic  Deformation  Theory 

The  following  algorithm  describes  an  iterative  procedure  by  which  the  anisotropic 
deformation  theory  could  be  implemented  into  a  code  like  ADINA.  There  are  many 
alternatives  that  are  simplier,  but  t he  method  described  below  appears  to  be  the  most 
straightforward. 

1)  Using  the  current  value  of  the  secant  modulus,  lCeP],  solve  for  the  new- 

total  displacements  in  the  next  load  or  time  step,  [Cc|1]  =  [Ce]  for 


2)  Using  a  Newton-Raphson  iterative  method  compute  the  new  effective 
stress,  a  ,  using  equations  ( C-2)-( C-7); 

3)  Calculate  the  new  estimate  of  the  components  of  the  plastic  strain  cP 
using  either  an  explicit  total  strain  or  an  explicit  total  work  method  (not 
described  herein); 

4)  Calculate  Sj,  rj,  Bjj,  \  and  the  total  plastic  strain  components  from  (C-tt) 
and  the  total  stresses  f mm  ( C-2a); 

5)  Repeat  steps  (2)-(ii)  until  convergence  is  achieved  for  a  .  and  . 
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